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Abstract 

We propose a phenomenological model of vortex dynamics in which the vor- 
tex is taken as a composite object made of two components: the vortex cur- 
rent which is massless and driven by the Lorentz force, and the vortex core 
which is massive and driven by the Magnus force. By combining the char- 
acteristics of the Gittleman-Rosenblum model (Phys. Rev. Lett. 16, 734 
(1966)) and Hsu's theory of vortex dynamics (Physica C 213,305 (1993)), 
the model provides a good description of recent far infrared measurements 
of the magneto-conductivity tensor of superconducting YBa2Cu307_5 films 
from 5 cm -1 to 200 cm -1 . 

PACS numbers : 74.25.Nf, 74.60.Ge, 73.50.Jt, 74.72.Bk 
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With the discovery of high-T c superconductors and the recognition of their potential for 
high field applications much attention has been devoted to the study of vortex dynamics. 
Nevertheless, the dissipation mechanism associated with vortex motion remains poorly un- 
derstood. From surface impedance measurements in the MHz and GHz ranges the vortex 
pinning and viscosity parameters have been reported based on the Gittleman-Rosenblum 
(G-R) model of vortex dynamics fl|-|3|. The observation of magneto-optical activity was re- 
ported in far infrared (FIR) experiments, which suggests a connection between vortex 
dynamics and cyclotron resonance ||. This observation greatly enriches the phenomenol- 
ogy of electrodynamics of vortices. Recently, quantized vortex core levels in YB^CuaOy^ 
(YBCO) have been probed directly using scanning tunneling microscopy f7j. A level spacing 
within the FIR frequency range was reported. Also evidence for vortex core excitations 
has been reported in FIR experiments ||, but the role of vortex core structure in vortex 
dynamics is not yet well understood. The use of terahertz spectroscopy P,|T0f and polarized 
FIR laser sources [II]] have recently filled in the gap of measurements between microwave 
and FIR frequencies. Therefore, a more satisfactory description of vortex dynamics over 
wide frequency range is needed, in order to describe these experimental observations. 

Microwave surface impedance studies [0-0,0 have demonstrated a loss band near zero 
frequency whose width is characterized by the depinning frequency Ud (= k/tj). There is 
remarkable consistency in the reported pinning force constant k for YBCO within a factor 
of two. However, the viscosity rj was found to vary much more widely. This leads to a 
wide range for ujd- From measurements of the conductivity up to 800GHz Ud ~ 300GHz 
(10 cm -1 ) has been reported for YBCO films ||. At FIR frequencies, optical activity has 
been observed in magneto-transmission measurements using broadband and laser sources 
throughout the 5 cm^ 1 < uj < 200 cm -1 frequency range PJ5|JTl[]. The conventional theories 
based on a Lorentz force acting on the vortices, e.g. the G-R model [|TI| and Coffey-Clem 
Model which are widely used to describe microwave surface impedance data, can not 
properly describe these FIR measurements since they do not produce optical activity ||. 
A model of vortex dynamics proposed by Hsu [14[] based on the microscopic BCS theory 
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in the clean limit gives a Magnus force acting on vortices and successfully describes the 
high frequency {u > 25 cm -1 ) optical activity in terms of hybridization of the cyclotron 
resonance and the vortex core resonance with the pinning resonance • The hybridized 



pinning-cyclotron resonance has been fully resolved in recent FIR laser measurements (u > 
5 cnr 1 ) |p,l|| . However, Hsu's theory fails to describe the low frequency absorption band JTTJ 
observed in the microwave studies |l|||[J. 

In this paper we present a new phenomenological model to describe the electrodynamics 
of type II superconductors. The two different inertial masses that have been discussed in 
vortex dynamics ]16| suggests a model consisting two coupled equations of motion which 
is intended to capture the qualities of G-R model at low frequencies and Hsu's theory 
at higher frequencies. We take the vortex to be a two component composite object (see 
Fig.|l]). One component is the vortex current outside the vortex core, which has a G-R type 
equation of motion, and the other component is the vortex core with an equation of motion 



based on Hsu's model [|TJ,|T^]. The resulting two coupled equations of motion produces a 
conductivity tensor with a London term and three Lorentzian oscillators : one low frequency 
oscillator similar to the G-R loss band and two finite frequency chiral oscillators as in Hsu's 
model. These three resonances are illustrated schematically in Fig.|2|. This composite vortex 
model provides a good description of the measured magneto-transmission and the magneto- 
conductivity tensor from 5 cm^ 1 < uj < 200 cm -1 |TT| . 

The motion of a vortex line in type-II superconductor under the influence of an alter- 
nating current is limited by frictional and pinning forces. The simplest description of linear 
vortex dynamics is ||12| 



J X O = f]V v + KT V (1) 

where J = nev s is the driving transport current density (v s is the superffuid velocity), v v 
(r v ) is the average vortex velocity (displacement), k is the pinning force constant, r\ is the 
viscosity and 0o is the flux quantum, hc/2e, which is expressed as a vector in the direction of 
the applied magnetic field (z). (A more precise interpretation of v v will be discussed later). 
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Eq. (|T[) does not contains an inertial term which corresponds to an assumption of massless 
dynamics of the vortex current pattern. This equation of motion produces a pole in the 
conductivity at zero frequency of width n/rj (Peak (2) in Fig.0) in addition to a London 
delta function reduced in strength from its zero field value. This model which is based on 
classical hydrodynamic analogies is valid only in those regions outside the core where the 
superfluid can be described by the order parameter in the time dependent Ginzberg-Landau 
(TDGL) theory [p!6|,|l7[ . Suhl [18| has considered the inertial mass of the vortex within 
the TDGL theory and concluded that it is order (A/e^) 2 smaller than the naive estimate 
M L ~ mniit; 2 . 

A Hall force term v v x z can be added to Eq. ([!]) leading to 

A (t> s — (5v v ) x z — rjv v — Kr v = (2) 

where A = nh/2 the strength of Lorentz force and f3 characterizes the relative strength of the 
Hall force to the Lorentz force. When f3 = 1 the first term corresponds to the Magnus force. 
A Hall term is required to describe the anomalous Hall effect observed in both conventional 
and high-T c superconductors in the mixed state where the sign of dc Hall effect changes if 
j3 changes its sign |T7]jT^ |. Otterlo et al. |Ll| have argued that, in a fermionic system, the 



strength of the Hall force depends on the electronic structure through the derivative of the 
density of states at the Fermi level. Their theory leads to a small (3 whose magnitude and 
sign depends on the details of the electronic band structure at the Fermi level. 

From an analysis of vortex motion within the BCS theory Hsu ]14|] proposed an equation 
of motion which has a kinetic term, vq corresponding to a non-negligeable vortex core 
mass, a Magnus force (v s — vc) X z and a 'gauge term' v s which couples the vortex motion 
to the external electric field. Ao and Thouless have shown that a Magnus force is a 
general property of vortex motion from a Berry phase argument. Hsu's equation of motion, 
normalized by the vortex core mass, is written as : 

vc = v s + Q (v s - v c ) x z - - v c - a 2 r c (3) 

T 



where hQ = A 2 ftp is the energy spacing of the quasi particle levels in the vortex core 
2TJ , vc [re] is the core velocity (displacement), - is the damping rate associated with 



the core motion and a is the bare pinning frequency. Due to the finite vortex core mass 
implicit to Eq. (|3|) two finite frequency poles of finite width are produced in the conductivity 
tensor in addition to a London term. One pole is a hybridized pinning-cyclotron resonance 
(Peak (1) in Fig.^j) and the other is a hybridized vortex core-pinning resonance (Peak (3) 
in Fig.|]). However, the vortex core resonance, which is quenched in the absence of pinning 
(leaving only the cyclotron resonance, in accord with Kohn's theorem ||), remains weak for 
physically reasonable pinning. 



The moving vortices produce a Josephson field p2] which couples to the superfluid by a 
modification of the London equation : 

eE 



m 



+ u c v L x z (4) 



where uj c = eH/ mc is the cyclotron resonance frequency and v £ is the average vortex velocity. 
Elementary arguments show that the Josephson field due to a moving vortex is proportional 
to the gradient of the local magnetic field distribution, E ~ V5(r). Since the curvature 
of the magnetic field is almost zero in the vortex core PB[, most of the contribution of the 



Josephson field should come from the current flow pattern outside the core, r ~ (London 
penetration length), where there is maximum curvature in the magnetic field distribution 
of the vortex. Consequently, we expect t£ to be the velocity of the vortex current pattern 
relative to the lattice. Therefore, 

vl = v v + v c . (5) 

where v v is interpreted as the velocity of the vortex current pattern relative to the core 
and Vc is the velocity of the core relative to the pinning center, as illustrated in Fig.|l]. In 
this picture the pinning force k r v acting on the vortex current pattern can be thought of 
as arising from the distortion of the vortex through its tendency to restore itself to the 
equilibrium configuration. Therefore, in this model k is an intrinsic property of the vortex. 
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As pointed out by Hsu [ll], there is an additional contribution to the current due to the 



moving vortex core which acts like a charged particle 

J = ne(Av s + (1 - A)v c ) (6) 

We determine A by requiring the response function to obey the energy conservation law. 
First we neglect the damping terms - vq and 77 v v and make the product E ■ J energy 
conserving. The only nonzero terms are the kinetic energy terms ^ (l^s) > J* (l^c 1 ) > j| \k v v) > 
and pinning potential terms ^ (| r c)5 ^ (l r ^)- Setting all other cross product terms to be 
zero, we obtain A = Q /(Q + u c ) which is similar to that obtained by Hsu |f[4}| . We also 
find that we have to add a term, — u c v v x z, to the right hand side of Eq.@. Combined 
with the 'gauge term' v s , v s — uj c v v x z = ^ + u c vc x z. The interpretation is that the 
total electric field acting on the vortex core is equal to the external electric field plus the 
Josephson field from the current flow pattern, but the Josephson field from the core itself is 
excluded. However, this is only a small correction since u c -C Q. 

From these equations, we can calculate the conductivity function. First from Eq.(f2|) and 
Eq.(|3p (in left circular polarization) we define 

g v (uj, H) = — = — — — — 7 

V s ipA + TJ + lK/UJ 

and 

v c -ico uj c g v (u,H) 
gc{u, H) = — = — ; — ——j 8 
v s —iuj + iil + ^ + ia' L /uj 

The conductivity function is then 

<r>, H) = l = ^ [Q + UJc 9 f c{uJ / H) \ /{Q ± Uc) , (9) 
E m -iuj - iuj c [gc{uJ,H) + g v {uj,H)\ 

from Eq.(^) and Eq.(^). o~ is obtained from a + through the time reversal symmetry relation 

a~ (u, H) = a + (—uj, H)* This relation suggests a canonical representation of a in which 

cr + is plotted for positive frequencies and a~ is plotted for negative frequencies. 

The conductivity function Eq.(g) is causal (i?e[cr] > 0), satisfies the oscillator strength 

sum rule, reduces to cyclotron resonance in the clean limit (in accord with Kohn's theorem 
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||) and further reduces to the London conductivity if the field is set to zero. With regard 
to the first two properties, we note that the conductivity tensor can be resolved into a 
finite sum of Lorentzian oscillators cr + (u, H) = ne 2 /mJ2l=o f q ~ u q) + ^q) ^ut with 

complex oscillator strengths, f q s, (except / the superfluid oscillator strength is real, see 
Eq.fllUD below.) and J2q=o fq = 1 as in G-R model and Hsu's theory. However, these 
complex oscillator strengths are unusual and may be an indication that the model needs 
further refinement. To obtain cyclotron resonance, we use the Bardeen-Stephen theory [24 



in which 77 = B&^o&n an d c n = r qp c 2 /47rA| where o~ n is the normal state conductivity, and 
r qp is the quasiparticle scattering rate. In the clean limit where r qp — > 00 and a 2 = 1/t = 
0, the conductivity reduces to —ne 2 /m /i(u — uj c ) ■ 

It is also interesting to consider the normal state limit of the model. In the flux flow 
regime where a 2 = but 1/t 7^ 0, and near T c the superfluid carrier is depleted so that A, 
Q — > (or <C 77, 1/r), we obtain a + (u,H) = ne 2 /m /[—i{u — u c ) + 1/t) which is the form 
of the metallic cyclotron resonance. However, a contribution from the quasiparticles should 
also be included, which would have the same form but proportional to n q = (n to tai — n) . 
Therefore, the model can have a smooth transition into the normal state. 

Recently, the FIR magneto-transmission coefficient, T ± (u,H) for YBCO thin films 
has been measured over a wide range of frequencies by the combination of broad- 
band Fourier Transform Spectroscopy (FTS) and CO2 pumped FIR laser measurements 



TT| . The magneto-conductivity tensor was obtained from the magneto-transmission 
spectra by a Kramers-Kronig analysis. Fig.^ shows the transmission amplitude ratio 
\T^(oj, H)/T(u, O)) 1 ^ 2 at H=9T and 4 K (Panel (a)) and the magneto-conductivity tensor 
Re[a + (uj,H)} (Panel (b)) in the canonical coordinate representation. A fit to the composite 
model (the dotted line) is also shown in Fig.|3] and is seen to be in very good agreement with 
the data except for some small structures above 50 cm -1 . These features are induced by the 
presence of 45° misaligned grains in these samples as discussed in Ref . The parameters of 
the best fit are a = 45 cm -1 , - = 58 cm -1 , Q = 49 cm -1 , uj c = 4.6 cm -1 P5fl , for the vortex 
core, and f3 = -0.086, k = 6xl0 5 N/m 2 , 77 = 5.6xl0~ 7 kg/m sec, for the vortex current, 
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assuming that A = nh/2 and ne 2 /m = c 2 /AitXq where Ao = 1850A ||11|| . By associating the 
Magnus force constant Myfl {My is the vortex mass) with nh/2 [pod , we can also estimate 
My pa 6 x 10 8 electron mass per cm. The corresponding estimate for the coherence length 
£v = 26A by taking My = mnirC, v |26|| . The fitting also gives Ud = k/t) ~ 10 cm -1 , very 
close to the value reported by Parks et al. |§. 

The resulting conductivity function is seen to be dominated by two major oscillators in 
addition to the London term (co> = 0, r = 0) |IT[] : a low frequency oscillator at ~3 cm -1 and 
width of 10 cm -1 , associated with the vortex currents; and a high frequency chiral oscillator 
at -24 cm -1 and width of 17 cm -1 in the hCP mode which is associated with the vortex 
core (the hybridized pinning-cyclotron resonance). Another oscillator has been identified at 
about 65 cm -1 in the eCP mode (the hybridized vortex core-pinning resonance), but it is 
very weak because of motional quenching || . 

The best fit also gives a negative (3 which may be related to the sign reversal of the dc Hall 
effect [|^.|l^.|l9|] . The presence of the Hall force (/3 ^ 0) in Eq.([|) shifts the low frequency 



oscillator to finite frequency uji ~ —j3ATi/r] (neglecting the contributions related to the 
vortex core, which is a fairly good approximation in this case.), where Ti ~ (k + Au c )/n 
is the width of this oscillator. (3 < results from uj\ > 0, which means that the low 



frequency oscillator is slightly electron-like (eCP). In the experimental paper |TTJ we also 
pointed out that the superconvergent sum rule, / °° In \T + (u, H)/T~(u, H)\ — = 0, which 
is a consequence of the supercurrent, is not satisfied unless the low frequency oscillator is 
shifted toward positive frequencies (eCP) to avoid the excessive low frequency weight in the 
optical activity. 

Finally we discuss fo the oscillator strength of the superfluid condensate within the 
composite vortex model. In the low frequency limit lim^^Q a + (u , H) = —fo(H) /(—icu), 
where 

1 + u) c (U/a 2 + A/k) 

Note that fo{H) = X 2 L /X 2 (H) which can be measured in microwave experiments and that 
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fo is independent of the damping parameters r\ and 1/r. In the G-R model, fo(H) = 
(1 + nhu c /2K,GR) 1 ■ Comparing Eq.(|T0|) with the G-R result we can define an effective 
11 k gr " in the low field limit where f (H) = 1 - cu c (tt/a 2 + A/k + Taking Myfi = A 

= nh/2 we find 1/"/cgk" = 1/k + 1/kc + 1/Afl where « is the force constant associated 
with the force between the vortex currents and the core and = Mya 2 is the force con- 
stant associated with the force between the core and the pinning center. In our model k 
is an intrinsic property of the vortex. A review of the published experiment data fl|-|3]|| 
indicates that the observed k is universal within a factor of two, which is difficult to un- 
derstand in terms of pinning to random defects in the lattice. The last term 1/Af2 is also 
intrinsic and comes from the reduction in the condensate density due to the vortex core 
volume (k/AQ = AirKm^y/nh 2 m 0.12). In general, we see that the small "force constant" 
dominates Jq. Therefore for strong pinning, such as occurs in YBCO films, the intrinsic pin- 
ning may determine u kgr" = k (in our case, k/kc ~ 0.13), which may explain the observed 
universal value reported in the literature. As gets small, for example in high quality 
single crystals, the extrinsic pinning Kc may dominate and u kgr" would become smaller 
and sample dependent. 

The authors acknowledged helpful discussions with T. C. Hsu, S. G. Kaplan, J. Orenstein, 
S. Wu, and V. Yakovenko. This work was supported in part by National Science Foundation 
under Grant No. DMR 9223217. 
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FIGURES 

FIG. 1. The vortex core (grey circle in the center) is pinned to the lattice defect (cross). The 
circulating vortex current pattern (concentric annular rings) is bound to the vortex core by the 
intrinsic restoring force (zigzag springs). 

FIG. 2. Three poles in the conductivity of the composite vortex model. Peak (1) and (3) come 
from the magneto-response of the vortex core. Peak (2) comes from the massless vortex current 
flow around the vortex core. The peak positions and oscillator strengths are the same as in Fig.^(b) 
but made much sharper for illustration purposes. hCP (eCP) stands for hole (electron) cyclotron 
resonance polarization mode. 

FIG. 3. (a) The transmission amplitude ratio \T^(u, H)/T(u, 0)|^ 2 as a function of frequency 
ijj for YBCO thin films at H=9T and 4 K. (b) The magneto-conductivity Re[a + (uj, H)] as a function 
of frequency u. The solid line represents the broadband data from 30 cm -1 to 200 cm -1 and the 
triangular points represent data from the laser source from 24.5 cm -1 down to 5.26 cm -1 . The 
dotted line is the fit to the composite vortex model. The dashed line in (b) is the residual metallic 
background Re[a(u>, 0)]. 
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